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Several methods of extrapolating the virial coefficients, including those proposed in this work, are discussed. 
The methods are demonstrated on predicting higher virial coefficients of one-component hard spheres. Esti- 
mated values of the eleventh to fifteenth virial coefficients are suggested. It has been speculated that the virial 
coefficients, B„, beyond B14 may decrease with increasing n, and may reach negative values at large n. The 
extrapolation techniques may be utilized in other fields of science where the art of extrapolation plays a role. 
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1. Introduction 

Virial coefficients play a key role in both molecular and phenomenological theory of fluids at low and 
medium densities. They are the coefficients in the density expansion of the equation of state expressed 
via the compressibility factor 

BP 00 

z=^— = i+£flip'- 1 , (i.i) 

P i=2 

where f} = l/ik^T) is the inverse temperature, P is the pressure, and p is the number density. Virial 
coefficients Bj are defined by exact formulae in terms of integrals whose integrands are products of 
Mayer functions (ljj. 

In this paper we concentrate on the simplest non-trivial model of fluid - the system of hard spheres. 
For hard spheres, i.e., hard body one-component fluids, the virial coefficients are numbers (they do not 
depend on temperature). For the system considered, the terms up to B4 are known analytically |2j,|3J]. It 
holds 

Bi = 1, (1.2) 

B 2 = 4, (1.3) 

B 3 = 10, (1.4) 

2707tt + 438 VI- 4131 arccos(l/3) 
B 4 = - = 18.364768... (1.5) 

707T 

The values are in the units of the packing fraction 77 = f c 3 p with a standing for the sphere diameter. 

The higher virial coefficients must be calculated numerically. Virial coefficients are sums of irre- 
ducible (cluster) integrals represented by diagrams |1]. It is a task for a high school student to determine 
the second virial coefficient for hard spheres. However, dimensionality of irreducible integrals increases 
rapidly with increasing order of the virial coefficient. In three dimensional space it is 3n - 6. For exam- 
ple the tenth virial coefficient contains 24-dimensional integrals (!) in the simplest case of spherically 
symmetric interactions between molecules. 
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Table 1. Number of unlabeled and labeled Mayer and Ree-Hoover diagrams, and dimensionality d of 
corresponding cluster integrals for hard spheres. 



n 


Mayer 


Ree-Hoover 


d 


unlabeled 


labeled 


unlabeled 


labeled 


2 


1 


1 


1 


1 


1 


3 


1 


1 


1 


1 


3 


4 


3 


10 


2 


4 


6 


5 


10 


238 


5 


68 


9 


6 


56 


11368 


23 


3053 


12 


7 


468 


1014888 


171 


297171 


15 


8 


7123 


166 537 616 


2 606 


56 671216 


18 


9 


194066 


50 680432112 


81564 


21286 987 064 


21 



Numbers of integrals (no general formula is known for the number of irreducible diagrams in de- 
pendence on ri) increase even much more rapidly than the dimensionality of the integrals, see table [l] 
Some of them are topologically equivalent which means that their integrands differ only in the num- 
bering of variables. In the table, the heading "labeled" denotes the total number of irreducible diagrams 
and "unlabeled" denotes the number of topologically different diagrams with numbered black points. To 
reduce these numbers, Ree and Hoover in their pioneering work |4] replaced Mayer diagrams with f- 
bonds by generalized Ree-Hoover diagrams with /-bonds and e-bonds using the identity f{r) + e(r) = 1. 
Besides the reduction of the number of topologically different diagrams, there is another advantage of 
the Ree-Hoover approach: computer codes using these diagrams are more efficient than the Mayer ap- 
proach. 

For higher virial coefficients, say n > 7, it is impossible within an average lifetime of an explorer 
to determine the diagrams and their weights using "pencil and paper" avoiding "human factor" errors. 
For example, for Bq there are 7 123 different Mayer diagrams which can be reduced to 2 606 Ree-Hover 
diagrams, and for each of them their weights should be determined. The analysis should be done auto- 
matically on a computer. Different algorithms utilizing symbolic algebra programming and automatic 
code generation have been proposed [5j-|8j]. 

Another problem is to evaluate the multi-fold integrals. This can be done by random-shooting Monte 
Carlo integration but it may require too much computer time. A more efficient way is to start from the 
so-called spanning diagrams - the integrals that can be calculated analytically. Configurations of the span- 
ning diagrams are sampled either using a standard Metropolis Monte Carlo method or (for the linear 
chains) by the so-called reptation |9]. The values of the diagrams of interest are then evaluated by Monte 
Carlo integration. 

Table 2. Summary of recommended virial coefficients B n for hard spheres. Values in parentheses are 
standard errors in the last significant digit. 



n 


B n 


5 


28.22445(10) 


6 


39.81545(34) 


7 


53.3418(15) 


8 


68.5394(87) 


9 


85.805(58) 


10 


105.8(4) 



The fifth virial coefficient for hard spheres was calculated by Rosenbluth and Rosenbluth (l^] and 
by Kratky jllUl3l . the sixth by Ree and Hoover Ql, the seventh by Ree and Hoover (3. by Ki m an d 
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Henderson (3» and by Janse van Rensburg and Torrie 1 17], the eighth by Janse van Rensburg 1 18] and 
by Vlasov et al. |5], the ninth by Labik et al. |6], and the tenth by Clisby and McCoy |7l|8j] . As a rule, while 
higher virial coefficients were just calculated, the lower ones were more precisely recalculated. The state- 
of-the-art values of the virial coefficients for hard spheres are shown in table[2]including their estimated 
uncertainties. 



2. Results and discussion 

2.1. Extrapolation of Virial Coefficients 

It is prohibitively difficult to calculate virial coefficients beyond Bio because the numbers of diagrams 
increase enormously. There are 2^1 diagrams with n field points that should be analyzed. It is 2 45 for 
n = 10, 2 55 for n = 11, 2 66 for n = 12, etc. The number of diagrams increases by more than three orders 
with n going from 9 to 10. It is thus practical to estimate higher virial coefficients for n > 10 rather than 
to try to calculate them. 

There are several approaches to extrapolation of virial coefficients based on an (unjustified) assump- 
tion that the higher coefficients depend on the lower ones. In the theory of fluids the most popular are 
Pade approximants, see 1 19] and references therein. They are based on the assumption that the com- 
pressibility factor may be expressed in the form 

Pniri) 

(2.1) 



where P n and Q m are polynomials of the order n and m, respectively, whose coefficients are combina- 
tions of the known virial coefficients. By expanding equation te.l) into the Taylor series, higher virial 
coefficients are obtained. 

Differential approximants (i(J are an extension of Pade approximants. The first order differential 
approximants have the form 

dZ 

PnW-r + Qm{r i )Z(ri) = R rn {T}), (2.2) 
ar\ 

where P n , Q m , and R m are polynomials; for P n (77) = 0, the Pade approximants are revealed. Differential 
approximants have not been used in the theories of fluids so far. 

These two approaches utilize only a part of available information, namely the lower virial coefficients. 
Another source of information are computer simulation data on the compressibility factors. They can be 
utilized by proposing an equation of state (EOS) some of its parameters being determined from the known 
virial coefficients and the remaining constants fitted to the simulation data. For example, Kolafa et al. EUl 
used an equation of state in the form proposed by Barboy and Gelbart [22|] 

" • 77 

Z=l+Y a t x l , x=—!—, (2.3) 

where lower ai are combinations of the known Bj and the higher ones are fitted to the computer simula- 
tion data. By expanding the equation in powers of 77, the higher virial coefficients may be estimated. 

This approach does not take into account the suggested convergence limit at high densities 1 23]. There- 
fore, we extended the form <2.3> by changing the pole 77 = 1 into 

n 

(2.4) 



*7pole - V 

where 77 po ie is a (nonlinear) adjustable parameter. This pronounced pole should not be confused with a 
weak nonanalyticity at the freezing point discussed below. 

All the traditional extrapolation methods suffer from an assumption that the values of lower virial 
coefficients are known precisely. This is, certainly, not true, see table [2] To be accurate, uncertainties in 
the lower virial coefficients (as well as in EOS simulation data) should be taken into account 1 24]. 
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We believe that the safest way of utilizing both the lower virial coefficients data and the EOS data is 
to minimize the objective function 



i=5 



(Bi-B exp \ 



+ E 



(Zj-Z' 



exp ,2 



a i 



(2.5) 



where n is a number of known virial coefficients, er,- are their standard deviations, k is the number of 
simulated state points, and a j are their standard deviations. Zj is the value of the compressibility factor 
at density r\ j given by a chosen compressibility factor in the form of equation 12.1) (Pade approximants), 
or in the form of equation d2.2t (differential approximants), or in the form of equation (2.3) (Barboy- 
Gelbart equation). 

As a measure of fit accuracy we use a standard deviation S 



(2.6) 



where p is a number of adjustable parameters. 



2.2. Barboy-Geibart approach 

Let us return to equations <2.3> and <2.4( . We allow for zero coefficients a,-s in the expansions. This im- 
proves precision without adding more adjustable parameters especially at higher densities. For instance, 
equation (10) of 12 ill , valid up to p — 0.98, uses powers x 1 to x 8 and x 12 . Therefore, the choice of powers 
serves as another (discrete) adjustable parameter. 

In order to elucidate the ability of expansions to predict unknown virial coefficients we repeated the 
calculations with forms d2.3> and <2.5) . using various sets of input data. That is, we "forget" that we know 
the virial coefficients up to B\q and use only a subset to n, 4 n *S 10. 

Moreover, we employ only subsets of available EOS data. For p = 0.94 the error in Bn ± 1 propagates 
to the error in Z + 0.00083 which is 8 times the standard deviation of the MD datum. The density range 
lower than p max < 0.94, would not contain enough information to allow for extrapolating the virial coeffi- 
cients. On the other hand, the equation of state is nonanalytical at the freezing point (density 0.947) |25], 
although the nonanalycity term is tiny and in moderately metastable region it is not detectable within 
available precision. Therefore, we choose p max = 0.98 as a compromise and calculate all fits with five 
EOS sets from low densities up to p max = 0.94, 0.95, 0.96, 0.97, 0.98. 



Table 3. Barboy-Geibart predictions of B n using the equation of state up to reduced densities 0.94 to 
0.98. Knowledge of all virial coefficients up to Bm was used, higher were predicted. The error bounds in 
parentheses include both systematic and statistical contributions. 



n 


M=4 


M=5 


M=6 


M=7 


M=8 f 


M=9 


M= 10 


5 


28.255(30) 














6 


39.67(7) 


39.716(23) 












7 


53.08(23) 


53.314(20) 


53.31(12) 










8 


69.0(2) 


69.36(9) 


68.6(2) 


68.66(5) 








9 


87.7(10) 


87.6(5) 


85.8(5) 


85.74(14) 


85.30(8) 






10 


109(3) 


106.8(8) 


105.9(18) 


104.9(2) 


105.1(5) 


106.3(1) 




11 


130(5) 


126.0(9) 


130.3(40) 


128(2) 


130.5(7) 


130.9(7) 


130.5(10) 


12 


150(3) 


142(2) 


160(5) 


158(7) 


164.2(9) 


156(3) 


156(3) 


13 


162(15) 


152(6) 


190(10) 


197(11) 


204(4) 


174(9) 


176(4) 


14 


160(40) 


159(20) 


230(40) 


230(15) 


232(26) 


175(10) 


184(8) 


15 


130(70) 


161(30) 


180(100) 


240(20) 


210(40) 


165(50) 


179(32) 



see figure[T] 
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Figure 1 . (Color online) Example of extrapolating the virial coefficients using formulas Z = 1 +L? j djX 1 + 
apxP based on B2 to Bq and EOS data up to p m ax- Virial coefficients B n , n > 8, are predicted, see the 
marked^ column of table [3] Black solid lines are the most precise fits (S < 1), blue dashed lines are less 
precise (S e [1, 1.2)), and red dotted lines the least precise (S e [1.2, 1.5)). The table is ordered according to 
the predicted S15 value and color-coded in the same way as the lines; some lines overlap. 



Table 4. The same as table[3]with expansion 12.4) . 



n 


M=4 


M=5 


M=6 


M=7 


M=8 f 


M=9 


M= 10 


5 


28.240(35) 














6 


39.73(11) 


39.72(2) 












7 


53.2(3) 


53.300(13) 


53.26(3) 










8 


69.1(3) 


69.24(7) 


68.50(3) 


68.67(3) 








9 


87.3(14) 


87.4(4) 


86.0(3) 


85.66(8) 


85.48(9) 






10 


107(4) 


107.0(7) 


106.4(7) 


104.80(16) 


105.3(3) 


106.15(6) 




11 


129(5) 


126.7(5) 


131.2(4) 


127.7(6) 


130.38(12) 


130.8(2) 


130.6(3) 


12 


148(9) 


144.7(8) 


162(3) 


157(4) 


162.2(24) 


158.1(14) 


158(2) 


13 


164(9) 


158(4) 


192(6) 


193(8) 


205(13) 


182(4) 


182(3) 


14 


174(35) 


165(12) 


207(16) 


227(11) 


234(35) 


192(7) 


192(4) 


15 


170(120) 


145(45) 


175(20) 


240(10) 


200(15) 


170(20) 


171(8) 



' see figure [T] 

The number of adjustable parameters is important in obtaining relevant predictions. We started with 
a small set and increased the number of adjustable parameters until the value of S dropped below 1.5. 
However, we never used more parameters if S < 1 was already reached to avoid a fitting noise. For 
each set of virial coefficients up to M, 4 M ^ 10, we analyzed several maximum powers in expansion 
d2.3> with all the above mentioned EOS data sets. The results were sorted into three categories: the best 
with S < 1, the good ones with S < 1.2, and the worst results with S < 1.5. Then, we plotted the curves 
of the predicted B n vs. n and graphically determined the median value and also tried to estimate the 
uncertainty from the data scattering. This error estimate is rather sensitive to the systematic errors which 
are larger than the statistical ones. We admit a subjective factor present in this procedure. One example 
is shown in figure [l] 

The results, based on more than a thousand of equations analyzed, are collected in tables [3] and |U 
It is seen that the procedure is capable of predicting the virial coefficients even though the errors are 
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sometimes underestimated. Of course, the virial coefficients obtained are less accurate than the directly 
obtained MC values. Only for Bio, the precision of the predicted value based on all previous coefficients 
approaches the precision of the MC result 0S- In addition, the accuracy of Bu and especially Bis de- 
creases. 

2.3. Pade approximants 

Simultaneous fitting of MD data and virial coefficients to rational functions (Pade approximants) is 
much more difficult because the equations are highly nonlinear. The approximants often give oscillating 
solutions from B12 or B13; we omit these solutions because we assume that the coefficients should behave 
"regularly". 

First we investigated the rational functions using all the known information, i.e., the exact values of 
B2 to B 4 , approximate values of B 5 toBio with standard deviations, and MD data to p max = 0.94, 0.95, 0.96, 
0.97, 0.98 (five tested variants). Approximants Pade [2,12], [8,5], [3,10], and [9,3] are the best (interestingly, 
[9,3] requires by 2 less parameters than [2,12]; [2,11] is worse, S = 1.4, pushing B14 to lower values) 
which with the five maximum densities gives 20 curves from which the predictions and their errors 
were determined. To test the predictive power of the procedure, we repeated the calculations without 
knowledge of B7 to Bio, see table [5] It is seen that the Pade approximants predict worse values of the 
lower virial coefficients than the Barboy and Gelbart approach. Nevertheless, the predictions of Bu to 
B15 using all the data are comparable. 

Table 5. Predictions using Pade approximants. See table[3]for the explanation of symbols. 



n 


M=6 


M= 10 


7 


53.28(7) 




8 


68.94(16) 




9 


84.7(14) 




10 


104.0(13) 




11 


130.6(4) 


130.5(10) 


12 


162(15) 


158(3) 


13 


195(30) 


181(4) 


14 


200(20) 


196(10) 


15 


150(60) 


198(24) 



2.4. Reliability of extrapolations 

While the statistical errors in the predicted higher virial coefficients can be easily evaluated, this 
cannot be said about the systematic errors introduced by a particular functional form. The only way to 
assess the systematic errors is to compare as many functional forms (Pade indices, powers in polynomials, 
etc.) as possible. Then, we subjectively estimate (after cancelling the evidently defective values) the "most 
probable" values of them. 

An example of a comparison of different methods of extrapolation is shown in table [6] for the four- 
teenth virial coefficient. In the table, ten particular extrapolation results are shown: three utilizing the 
Barboy-Gelbart equation 12. 3> (with and without adjustable virial coefficients), Pade approximants, and 
a differential approximant, equation <2.2> . 

The final recommended values of higher virial coefficients including the estimated uncertainties are 
given in table[7] They are based on critical examination of several methods and many equations of state. 
In the table, our values are compared with recent literature estimates obtained independently using dif- 
ferent approaches by different authors 171 18l l26l429i] . Except B15 they agree with ours within estimated 
errors. However, a general trend is that we predict Bu to B13 higher whereas B14 and B15 lower. 

In the table our estimated values of B n are shown up to n — 15. We calculated higher virial coefficients 
as well. However, their values exhibited a large scatter and no final estimates were thus made. Anyway, 
the literature estimates for B\e are given in the table. 
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Table 6. The extrapolated fourteenth virial coefficient using all the known virial coefficients and densities 
up to p = 0.94. BG+1 denotes the Barboy-Gelbart equation with one additional term x 10 , BG+11 with 
x 10 and x , etc., and no optimization of the virial coefficients values, 0- in front of symbols means that 
the full optimization, {2.5\ , was performed. Pade and O-Pade denote the nonoptimized and optimized 
Pade approximants, equation (2. It . O-DA is the optimized differential approximant. 



Method 


B14 


S 


BG+1 


-270 


234 


BG+11 


55 


19 


BG+111 


198 


1.54 


O-BG+1 


215 


1.59 


O-BG+11 


193 


1.45 


O-BG+111 


186 


0.97 


Pade [5,4] 


204 


8.3 


Pade [2,7] 


198 


11.3 


O-Pade [5,4] 


207 


6.9 


O-Pade [2,7] 


206 


6.7 


O-Pade [9,1] 


204 


6.5 


O-Pade [7,3] 


201 


6.3 


O-DA 


200 


6.3 



Table 7. Extrapolated virial coefficients. 



n 


this work 


ref. [7, 8] 


ref. [26] 


ref. [27] 


ref. [28] 


ref. [29] 


11 


130.6(10) 


128 


128 


129 


126 


128 


12 


137(3) 


153 


153 


154 


149 


153 


13 


183(10) 


181 


181 


182 


174 


182 


14 


190(20) 


215 


213 


211 


202 


214 


15 


180(30) 


247 


248 


236 


231 


247 


16 




279 


288 


254 


262 


279 



We cannot assert that really B\s < Bu; in fact, oscillating behavior is typical of polynomial approxi- 
mants and these results are consistent with a monotonous increase, at least for not too large n. Ultimately 
the sequence of B n cannot be monotonous because the equation of state is not analytical at the freezing 
point and consequently its radius of convergence cannot exceed the freezing density (we believe that it 
equals the freezing density). 

3. Concluding remarks 

In this work we have used several extrapolation methods to estimate the unknown virial coefficients 
of hard spheres from B\\ to B\$ and compared them with recent literature estimates. We believe that the 
presented results could be used in testing theoretical approaches of the thermodynamics of hard spheres. 

The problem of reliability of any extrapolation result rests in the estimate of its uncertainty. This is, 
in principle, subjective depending on an author's view (either optimistic or pessimistic). Here we have 
tried neither to be too optimistic nor too pessimistic. We have used comparisons of extrapolated virial 
coefficients obtained using as many different approaches as possible. 

The techniques used in this work can be readily extended to virial coefficients of hard body fluids 
such as spherocylinders, diatomics and multi-atomics, mixtures of hard spheres and other hard body 
fluids and their mixtures. What one needs is i) reliable values of lower virial coefficients and estimates of 
their uncertainties, and ii) precise equation-of-state computer simulation data. An extrapolation of virial 
coefficients of non-hard-body fluids is less straightforward due to their dependence on temperature. 



23004-7 



M. Oncak etal. 



References 

1. Hill T.L., Statistical Mechanics, McGraw-Hill, New York, 1956. 

2. Nijboer B.R.A., Van Hove L., Phys. Rev., 1952, 85, 777; doi 10.1103/Phys Rev.85.777| 

3. Lyberg I., e-print cond-mat/0410080. 

4. Ree F.H., Hoover W.G., J. Chem. Phys., 1964, 40, 939; doi |10.1063/1.1725286| 

5. Vlasov A.Y., You X.-M., Masters A.J., Mol. Phys., 2002, 100, 3313; doi 10.1080/00268970210153754 

6. Labik S., Kolafa J., Malijevsky A., Phys. Rev. E, 2005, 71, 021105; doi 10.1103/PhysRevE.71.021105 

7. Clisby N., McCoy B.M., J. Stat. Phys., 2006, 122, 15; doi |10.1007 /sl0955-005-8080^0l 

8. Clisby N., McCoy B.M., Pramana-J. Phys., 2005, 64, 775; doi 10.1007/BF02704582 

9. Allen M.P., Tildesley D.J., Computer Simulation of Liquids. Clarendon Press, Oxford, 1987. 

10. Rosenbluth M.N., Rosenbluth A.W., J. Chem. Phys., 1955, 23, 356; doi 10.1063/ 117419671 

11. Kratky K.W., Physica A, 1976, 85, 607; doi 10.1016/0378-4371(76)90029-7 

12. Kratky K.W., Physica A, 1977, 87, 584; doi 10.1016/0378-4371(77)90051-6 

13. Kratky K.W., J. Stat. Phys., 1982, 27, 533; doi 10.1007/BF01011091 

14. Kratky K.W., J. Stat. Phys., 1982 29, 129; doi 10.1007/BF01008253 

15. Ree F.H., Hoover W.G., J. Chem. Phys., 1967, 46, 4181; doi 10.1063/1.1840521 

16. Kim S., Henderson D., Phys. Lett. A, 1968, 27, 378; doi 10.1016/0375-9601(68)91066-9 

17. van Rensburg E.J., Torrie G.M., J. Phys. A: Math. Gen., 1993, 26, 943; doi: 10.1088/0305-4470/26/4/022 

18. van Rensburg E.J., J. Phys. A: Math. Gen., 1993, 26, 4805; doi: 10.1088/0305-4470/26/19/014 

19. Sanchez I.C., J. Chem. Phys., 1994, 101, 7003; doi |10.1063 /1.468456 

20. Guttmann A.J. - In: Phase Transitions and Critical Phenomena, Vol. 3, chapter 1, eds. C. Domb and J. Lebowitz. 
Academic Press, 1989. 

21. Kolafa J., Labik S., Malijevsky A., Phys. Chem. Chem. Phys., 2004, 6, 2335; doi |10.1039/b402792b| 

22. Barboy B., Gelbart W.M., J. Chem. Phys., 1979, 71, 3053; doi 10.1063/1.438711 

23. Yelash L.V., Kraska T., Phys. Chem. Chem. Phys., 2001, 3, 3114; doi 10.1039/bl02972j 

24. Erpenbeck J.J., Wood W.W., J. Stat. Phys., 1984, 35, 321; doi |10.100 7/BF01014387l 

25. Kolafa J., Phys. Chem. Chem. Phys., 2006, 8, 464; doi 10.1039/b511 999el 

26. Tian J.X., Gui Y.X., Mulerlo A., Phys. Chem. Chem. Phys., 2009, 11, 11213; doi[l0Tl039/b915002al 

27. Bannerman M.N., Lue L., Woodcock L.V, J. Chem. Phys., 2010, 132, 084507; doi 10.1063/1.3328823 

28. Hu J.W., Yu Y.X., Phys. Chem. Chem. Phys., 2009, 11, 9382; doi 10.1039/b911901al 

29. Tian J.X., Gui Y.X., Mulerlo A., J. Phys. Chem. B, 2010, 114, 13399; doi 10.1021/jpl06502x 



flo excTpano/isiLjiT Bipia/ibHi/ix i<oe<t>ii^ieHTiB TBepflux c<|>ep 

M. Cd-waKR A. Ma/iieBCKP, M. Ko/ia^aPl C. /la6iP 

iHCTHTyT xiMinHoT TexHO/iorii, npara, Biflfli/i cf>i3i/iHHOi xiMii, 166 28 npara 6, MecbKa Pecny6yiiKa 

Biflfliyi <j)i3HKH, (jjaKy/ibTeT npwpoflHW4wx HayK, yHiBepci/ueT m. OcrpaBa, 701 03 OcrpaBa 1, HecbKa Pecny6/iiKa 

06roBopK)K)TbCfl fleKi/ibKa eKCTano/iflUMHwx MeTOfliB A/in Bipia^bHux Koecpii4ieHTiB, bk/ikjhho i3 3anponoHO- 

BaHUMH B UM po60Ti. riOKa3aHO MO>K/1MBOCTi MeTOfliB nepefl6a4aTH BWLUi Bipia^bHi KOetjjiL^iEHTH OflHOKOMnO- 

HeHThtoi cucTeiviw TBepflHx cefiep. 3anponoHOBaHO oi4iHKn pj\n Bipia^bHwx Koe(()iL4ieHTiB Bifl oflUHafl^Toro pp 
n'flTHafl^Toro. nepeA6a4aeTbC5i, mo Bipia/ibHi Koe((iii4ieHTM B n Bnu(oro nopaflKy Hix B14, Moxyrb 3MeHLuyBa- 
™cfl 3 pocTOM n i MOxyTb floc^raTH HeraTHBHUx 3Ha4eHb pj\n Be/11/iKwx n. \\\ eKCTpano/ini4itiHi ue\opy\ MOxyTb 
6yTH BUKopucTaHi b iHLUnx 06/iacTflx, p,e eKCTpano;ifli4ia e Ba>icfinBOK). 

ICnioMOBi c/ioBa: iBepp) ccpepn, BipianbHi Koetpiiiiei-m/i, eKcrpano/ifiLfiMHi Merofli/i 
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